arXiv: 1506.00341 vl [cond-mat.quant-gas] 1 Jun 2015 


Quantum Crystals and Laughlin Droplets of Cavity Rydberg Polaritons 


Ariel Sommer,^ Hans Peter Biichler,^ and Jonathan Simon^ 

^Department of Physics and James Franck Institute, University of Chicago, Chicago, IL 
^Institute for Theoretical Physics III, University of Stuttgart, Cermany 
(Dated: June 2, 2015) 

Synthetic quantum materials offer an exciting opportunity to explore quantum many-body physics 
and novel states of matter under controlled conditions. In particular, they provide an avenue to 
exchange the short length scales and large energy scales of the solid state for an engineered system 
with better control over the system Hamiltonian, more accurate state preparation, and higher fidelity 
state readout. Here we propose a unique platform to study quantum phases of strongly interacting 
photons. We introduce ideas for controlling the dynamics of individual photons by manipulating 
the geometry of a multimode optical cavity, and combine them with recently established techniques 
to mediate strong interactions between photons using Rydberg atoms. We demonstrate that this 
approach gives rise to crystalline- and fractional quantum Hall- states of light, opening the door to 
studies of strongly correlated quantum many-body physics in a photonic material. 


Each of the synthetic quantum material platforms 
currently under exploration provides an exquisite win¬ 
dow into the physics of condensed matter: ultracold 
atoms in optical lattices mm benefit from strong in¬ 
teractions, extraordinary coherence, high-fidelity read¬ 
out mg, and excellent optical control over lattice ge¬ 
ometry [3, enabling studies of the Bose- and Fermi- 
Hubbard models mm, along with precision measure¬ 
ments of bulk physics H 11 [g. Trapped ions have 
been employed for small- to medium- scale studies of 
quantum magnetism US E], and classical crystalliza¬ 
tion [laiis]. Recently, arrays of coupled microwave res¬ 
onators have been explored as a promising platform for 
high-fidelity quantum simulation of dynamical equilib¬ 
rium lattice physics [T4HT8] . Reaching a regime of strong, 
long range interactions, synthetic magnetic fields, and 
real-time particle injection would open unique routes to 
probe the interplay of emergent crystallinity, topology, 
dynamics and dissipative state preparation in a regime 
inaccessible in other materials. This work describes a 
route to these objectives in a two-dimensional platform 
combining multi-mode optical cavities to control single¬ 
particle dynamics, with Rydberg slow light polaritons to 
induce strong interactions. 

The combination of electromagnetically induced trans¬ 
parency (EIT) with the strong interaction between Ry¬ 
dberg atoms has recently emerged as a tool to induce 
strong interactions between individual photons [I^27j . 
The appearance of a bound state of photons [ 55 ] and the 
realization of single photon transistors [551 HO] demon¬ 
strate the power of this approach, while theoretical pro¬ 
posals have explored the extensions to two-qubit pho¬ 
tonic gates m isiHsa. In the one-dimensional free- 
space setups realized so far, the optical depth limits the 
number of interactions per photon. Recently, a single¬ 
mode optical cavity was employed to enhance the photon- 
photon collision probability per photon lifetime, lead¬ 
ing to a dispersive nonlinearity [551 |5g. Meanwhile, 
theoretical proposals have begun to explore arrays of 
coupled single-mode cavities to realize photonic lattice 
physics HU Hg. On the other hand, multimode opti¬ 


cal cavities have been proposed as a platform to simulate 
glassy physics Hi HZ], and used to observe Bose-Einstein 
condensation of photons thermalized with dye |38j or hy¬ 
bridized with excitons [39ti4T] . These systems operate in 
an open limit where mean-field dynamics of weakly in¬ 
teracting photons can be observed in real-time and par¬ 
ticles are injected into the system as desired. It is sensi¬ 
ble to consider marrying multimode optical cavities with 
the strong photon-photon interactions accessible through 
Rydberg EIT. We explore this possibility, and develop a 
simple, physical framework for understanding the result¬ 
ing quantum many-body system. 



FIG. 1. Rydberg polaritons in a multimode optical 
cavity, a, A quasi-two-dimensional atomic gas of ground- 
state alkali atoms is held in the waist of an optical cavity 
formed by a set of mirrors. Photons enter the cavity through 
a weak driving laser (red) and interact through hybridization 
with atomic Rydberg excitations, induced by a strong control 
laser (blue). Information about the state is read out by detect¬ 
ing the transmitted photons, b, Relevant atomic level struc¬ 
ture giving rise to Rydberg polaritons through electromagnet¬ 
ically induced transparency (EIT). Cavity photons couple the 
atomic ground state \g) to an intermediate state |e) with col¬ 
lective Rabi frequency G, while a strong control laser couples 
|e) to a Rydberg level |r) with Rabi frequency Cl. The con¬ 
trol laser is detuned by A from single photon resonance, and 
together with the cavity frequency gives a mode-dependent 
two-photon detuning of 5 from EIT resonance. 
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Our proposed approach employs photons in a family 
of near-degenerate resonator modes to mimic the physics 
of a two dimensional gas of massive particles in a trap. 
The modes must be nearly degenerate so that photons 
can be coupled between them via Rydberg-mediated in¬ 
teractions, giving rise to photon-photon collisions. In 
practice, such a setup would consist of a high-finesse op¬ 
tical cavity to engineer the photonic modes, along with 
a gas of laser-cooled atoms within the cavity to mediate 
photon-photon interactions (Fig. [^. In what follows we 
provide a formalism that describes the dynamics of the 
photons as massive trapped particles in the presence of 
synthetic magnetic fields; introduce coupling to the Ry¬ 
dberg EIT medium, and compute a renormalized photon 
mass and interparticle potential; and perform numeri¬ 
cal experiments demonstrating that few-body phenom¬ 
ena such as crystallization and Laughlin droplet forma¬ 
tion are directly observable in such a system. 


I. CAVITY PHOTONS AS PARTICLES 

It is a remarkable property of nearly degenerate optical 
cavities that the behaviour of a photon in the transverse 
plane is well described as a massive particle in an exter¬ 
nal potential in two dimensions. This behaviour can be 
understood in both geometric and wave optics pictures. 
We begin with geometric optics to provide intuition. 

In a degenerate cavity, precise tuning of the geome¬ 
try causes light rays to retrace the same paths repeat¬ 
edly, resulting in a fixed set of intersection points with a 
transverse plane. Tuning slightly away from degeneracy 
leads to imperfect repetition of the ray paths, causing 
the intersection points to process and trace out the path 
of a particle moving in two dimensions. In the case of 
a near-planar Fabry-Perot cavity with spherical mirrors, 
where the radius of curvature of the mirrors greatly ex¬ 
ceeds the distance between the mirrors, the particle ex¬ 
ecutes harmonic oscillation, as show in Fig. [^. The 
transverse oscillation frequency a;_L depends on the cav¬ 
ity geometry, and is independent of the light wavelength. 
Near degeneracy, wj, becomes much smaller than the fre¬ 
quency of round trips, allowing one to coarse-grain over 
the longitudinal motion and consider purely transverse 
two-dimensional motion. 

In the wave optics picture, a photon lives in modes 
with three quantum numbers: one longitudinal and two 
transverse. For cavities that are short compared with 
their mirror radii of curvature, the transverse mode spac¬ 
ing is much smaller than the longitudinal mode spacing, 
so coarse-graining over the optical round-trip amounts 
to considering only a fixed longitudinal quantum num¬ 
ber. In particular, two-mirror Fabry-Perot cavities ex¬ 
hibit Hermite-Gauss (HG) eigenmodes [ 32 ] j with frequen¬ 
cies uJmn = uj±(m+n+l)+const for the transverse modes 
HGmn, illustrated in Fig. In the transverse plane of 
the cavity waist, the HG modes have the same form as 
the eigenfunctions of the two-dimensional harmonic os- 



FIG. 2. Photons in near-degenerate cavities as parti¬ 
cles in two dimensions, a, Light rays (blue) intersect the 
transverse plane of the cavity at slightly different locations 
after each round trip, owing to the curvature of the mirrors. 
The intersection points trace out the motion of an effective 
particle undergoing harmonic oscillation. The cavity eigen¬ 
modes in the wave picture correspond to the wavefunction 
of the particle, b, Mode spectrum of a near-planar cavity 
within a manifold of transverse modes. The splitting between 
levels is small compared to the cavity free spectral range (not 
shown), and approaches zero in the planar limit, c, Schematic 
of a helical running-wave cavity that induces an image rota¬ 
tion by an angle 4>- d. Frequency spectrum in a single lon¬ 
gitudinal mode of a helical, showing the shift of the mode 
frequencies proportional to the angular momentum 1. Modes 
with negative I are shifted upward in frequency, while modes 
with positive I are shifted downward towards degeneracy with 
the TEMoo mode, forming a Landau level. 


dilator, with an oscillator length of wq/V^, where wq 
is the cavity waist size (1/e^ intensity radius). This 
energy- and mode spectrum corresponds to that of the 
quantum harmonic oscillator, thus the photons in the 
cavity may be viewed as 2D particles near the quantum 
ground state of a harmonic trap. The photon “mass” 
then arises from the analogy to zero-point motion, corre¬ 
sponding to the 1/e^ intensity radius of the lowest cavity 
mode Wq = ^y/L{2R — L), where L is the cavity length, 
R is the radius of curvature of the cavity mirrors, and A 
is the wavelength of the light. The photon mass is then 
ruph = 2H/{wqUJj_) for HG modes. In the special case of 
a near-planar cavity, the mass reduces to the relativistic 
expression ruph = huj±fc^, with c the speed of light [38] . 

More exotic cavity geometries give rise to more com¬ 
plex transverse dynamics of the photons in the focal 
plane. Of particular interest are geometries that pro¬ 
duce an image rotation on each round-trip, arising from 
a running-wave geometry with either intra-cavity dove 
prisms, or non-planar geometry (Fig. [2:). Neglect¬ 
ing astigmatism, such helical cavities exhibit Laguerre- 
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Gauss (LG) eigenmodes that carry orbital angular mo¬ 
mentum hi in the transverse plane, where I is an inte¬ 
ger. The modes are shifted in frequency by with 

Wrot = Lr4’|c^ where </> is the round-trip rotation angle 
and Lr is the round-trip distance. Tuning Wrot therefore 
brings modes with different I into degeneracy, as illus¬ 
trated in Figure [^. The LG modes can be indexed by 
positive integers m and n counting units of positive and 
negative angular momentum |43j so that I = m — n. At 
degeneracy, the frequency spectrum is independent of one 
of the indices (m). The degenerate manifolds correspond 
to Landau levels, where the lowest Landau level has n = 0 
and consists of the transverse modes 


'^mo{p) 





( 1 ) 


Here p = (x, y) is the transverse position and z = x + iy. 
The magnetic length Is = wo/2 sets the product of the 
effective charge and effective magnetic field to AH/wq, 
giving 4 flux quanta per mode area (ttwq). The cy¬ 
clotron frequency Wc equals the frequency spacing be¬ 
tween Landau levels, determined by the cavity geome¬ 
try, and determines the mass ruph = Ah/( wqUJc)- Tuning 
slightly away from degeneracy, the frequency spectrum 
becomes ujmn = muj'^/ujc + nuic, and induces, in addi¬ 
tion to the magnetic field, a harmonic potential with 
frequency a;_L that vanishes at degeneracy. The com¬ 
plete coarse-grained Hamiltonian governing the photon 
dynamics is thus: 

hph = {-ihdp - A)^ -b ^WphW^p^ (2) 

zruph I 

with A = AB{—y,x,0). 

In more general degenerate cavities, light rays only re¬ 
trace their paths after s > 1 round trips. Goarse-graining 
must then also incorporate multiple round trips, result¬ 
ing in near-degenerate manifolds with mixed longitudinal 
quantum numbers. The transverse modes in such near¬ 
degenerate manifolds span only a subset of the complete 
Hilbert space, and the resulting photon dynamics exhibit 
additional symmetries in real- or phase- space. In the 
case of helical cavities, degeneracies with s > 1 corre¬ 
spond to particle motion on the surface of a cone |44j. 
The Methods section provides additional details on the 
particle description of cavity photons for general degen¬ 
eracies. 


II. COUPLING TO A RYDBERG EIT MEDIUM 


to a very large principal quantum number n ^ 100 inter¬ 
act very strongly [45]. An emerging technology laiEi] 
hybridizes photons with Rydberg excitations of ground- 
state atoms to produce photonic quasiparticles whose in¬ 
teractions come from their Rydberg part, and motion 
from their photonic part. A second (strong) laser beam 
helps with this hybridization, as the ground-to-Rydberg 
oscillator strength is very small, making direct absorption 
to the Rydberg state very unlikely. The probe photons 
are absorbed with high probability on a low-lying, nearly 
closed atomic transition, and from which the atoms are 
rapidly excited to the Rydberg state by the strong beam. 
In steady state, these couplings lead to EIT, in which a 
probe photon propagates as a dark polariton [46]. The 
Rydberg component of the dark polariton then leads to 
strong polariton-polariton interactions. 

We propose using the Rydberg EIT technique de¬ 
scribed above to hybridize cavity photons with Rydberg 
excitations of an atomic gas held in the cavity (Fig. 

The atomic gas-we consider ®^Rb-is confined to 
the transverse plane of the cavity waist in a layer that 
is thin compared to the Rayleigh range zr of the cavity 
waist so that longitudinal diffraction of the photons may 
be neglected; this will be important to ensure exclusively 
real-space two-body interactions, as photon diffraction 
corresponds to a fractional fourier transform [?7]. We 
further require that the sample be sufficiently optically 
thin as to not mix longitudinal modes of the optical res¬ 
onator. 

A strong control laser couples the first electronic ex¬ 
cited state of the atoms |e) to a Rydberg level |r) with 
Rabi frequency fl and detuning A, as shown in Fig. 
Cavity photons are tuned near EIT resonance, allowing 
them to propagate in the medium as Rydberg polaritons. 
In the transverse plane, the atomic medium has a uniform 
density and extends out to a radius that greatly exceeds 
the cavity waist wq. Consequently, the atomic medium 
does not mix different transverse modes of the cavity, 
except through the induced photon-photon interaction. 
Individual polaritons therefore have the same qualitative 
properties as photons in the bare cavity, with rescaled 
parameters. The Hamiltonian describing slow light Ry¬ 
dberg polaritons inside a multimode cavity is given by 
(see Appendix for details of the derivation): 


iL = iLph + iLeit + Hint (3) 


in which we include non-Hermitian terms to account for 
dissipation. The first term contains the photon dynamics 
introduced in the previous section. 


Photons do not interact directly with one another, 
and therefore a photonic material requires a nonlinear 
medium to mediate photon-photon interactions. Most 
optical nonlinearities are weak at the few-photon level, 
making them unsuited to creating strongly-correlated 
photonic materials. On the other hand, atoms excited 



where '0a is the two-dimensional field operator for pho¬ 
tons in a nearly degenerate manifold. The second term 
in (§ contains the coupling of cavity photons to collec¬ 
tive atomic excitations with Rabi frequency G, and the 
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coupling to the Rydberg level, 

/ -iK G 0 \ /ijjaX 
G 2A n V'e (5) 

y 0 n -17 y yV'T-/ 

where and tpr are the two-dimensional field operators 
for collective excitations to the atomic excited state and 
Rydberg state, respectively. Here k is the cavity loss 
rate, A = A — ir/2 is the complex single-photon detun¬ 
ing, where F is the atomic excited state decay rate, and 
7 is the Rydberg level decay rate. The third term in 
^ describes the interactions between pairs of Rydberg 
atoms through a potential C/(r), 

Hint = ^ J drdr'U{r - r')(j)l{r)(j)l{r')(j)r{r)(l)r{r') ( 6 ) 

where is the field operator for local excitations 

of the Rydberg level in three dimensions. We con¬ 
sider Rydberg S levels giving a van der Waals potential 
U{r) = —CqIv^. At high polariton density, additional 
three-body interactions can appear [48l|49j. Finally, pho¬ 
tons enter the cavity through driving by a weak laser at 
a frequency Sl relative to FIT resonance, which spectro¬ 
scopically probes the Hamiltonian ([^ and drives it to a 
dynamical steady state. 

The low-energy excitations of the system consist of 
dark polariton quasiparticles gSlIiS]- A dark polariton 
has probability cos^ 6* to be a photon and a probability 
sin^ 0 to be a Rydberg excitation, where 9 = tan“^(G/H) 
gives the dark state rotation angle. The dark-polariton- 
projected single-particle Hamiltonian is given by a rescal¬ 
ing of the photon Hamiltonian (see Appendix for details 
of the derivation), 

^dark — ^ph COS 9 

= iy-ihdp - Af + ^mDUjl_DP^ ( 7 ) 

where the polariton oscillation frequency is w_ld = 
w_L cos^ 9 and the mass \s mo = niph/ cos^ 9. 

The Hamiltonian 0 has discrete eigenstates that each 
correspond to a dark polariton in a particular cavity 
mode. The eigenfrequency of a dark polariton, measured 
relative to the EIT resonance frequency for a probe pho¬ 
ton, is then given by S cos^ 9, where S is the detuning of 
the cavity mode from EIT resonance. Likewise, the loss 
rate of a dark polariton state is given to leading order by 
K cos^ 0-1-7 sin^ 9, corresponding to out-coupling through 
the cavity mirrors and loss through decay of the Rydberg 
state. For larger detunings, breaking of the EIT condi¬ 
tion leads to additional loss proportional to S'^. Figure 
ir -b shows the frequencies of the single-polariton eigen¬ 
states of the system, including the dark polariton state 
as well as four higher-frequency states. 

The Rydberg-Rydberg interactions detune the Ryd¬ 
berg level at short range, leading to a screening of the 
divergence of U{r) at r = 0. The screening results 
from mixing of the higher-frequency single-particle states 




r (urn) 


r (pm) 


FIG. 3. Dispersion and interaction of cavity Rydberg 
polaritons. a, Eigenfrequencies of the single-polariton states 
coupled to a cavity mode with detuning 5 from EIT resonance. 
Red: dark state, black: bright states, gray: matter-only 
states, independent of 5, in spatial modes not coupled to the 
cavity, b. Magnified view of the dark polariton frequency Sd- 
The lightly shaded region shows the dark polariton linewidth 
(imaginary part of the frequency). The dark shaded region 
shows the contribution of the photon out-coupling, c. Ef¬ 
fective potential Ud between two dark polaritons versus dis¬ 
tance r. The short-range limit is given by the parameter 
— 1/x defined in the text, while the long-range behaviour is 
proportional to 1/r® for van der Waals interactions. Solid: 
real part, dashed: imaginary part, d, Rydberg-Rydberg com¬ 
ponent (firr of the wavefunction for two dark polaritons in 
a single-mode cavity. The extent of the wavefunction is de¬ 
termined by the cavity waist wo, while the wavefunction is 
suppressed at short range r < ^, with ^ the blockade radius. 
The parameters are G = 27r • 27 MHz, re = 27r • 1.0 MHz, 
H = 27r • 5 MHz, n = 10^^ cm“®, wo = 30/rm, Lr = 6.8 cm, 
and Ce = —27r • 4243 GHz pm® for the SOS Rydberg level of 
®"Rb. 


with the dark state in the two-body wavefunction. At 
low energies, the screened potential becomes C/®®(r) = 
-Ge/iCeX + as in Ref. [57], where x is given by 

^ 2(H/2)2 " 2A 

and provides a characteristic length scale ^ 
for the blockade of two Rydberg excitations. The inter¬ 
action potential Ud(,p) between dark polaritons is then 
determined by restricting to two dimensions and pro¬ 
jecting onto the dark polariton basis, as we show in 
the Appendix. In the limit where the atomic medium 
is thin compared to the blockade radius, it reduces to 
Ud ~ <aC/®®(/9), where a = sin"* 0 is the probability for 
two dark polaritons to both be Rydberg excitations. Fig¬ 
ure 0 shows the interaction potential for a realistic set of 
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FIG. 4. Few-body emergence of crystalline structures. The restriction to even harmonic oscillator wavefunctions for a 
cavity with s = 2 leads to unique crystals where individual photons occupy a superposition of two lattice sites, a, Illustration 
of the photon probability for four photons in a cavity with s = 2. The red dots mark possible positions of the photons in 
the cavity in a single shot, b and c show numerical calculations for two photons in a cavity with one transverse dimension, 
implemented using TEMmO modes, b, Probability to find a photon versus transverse position in the cavity. One photon (red) 
is predominantly sitting in the center of the trap, while the second photon (blue) can be found with equal probability either to 
the left or to the right. The dashed line illustrates the lowest harmonic oscillator wavefunction, while the green line indicates 
the effective interaction potential, c, Fluid-to-crystal transition of two photons with increasing interaction strength. The color 
scale shows the emitted light intensity distribution in arbitrary units. 


parameters, while Fig. P shows the two-Rydberg part 
of the wavefunction in a two-polariton state. 

The length and time scales governing the two- 
dimensional gas of Rydberg polaritions allow direct mea¬ 
surements on the system using optical techniques. Con¬ 
ventional optics can resolve the blockade radius of typi¬ 
cally lO^m, which sets the smallest distance between po- 
laritons in a stable configuration. Likewise, existing elec¬ 
tronics can resolve the timescale set by the FIT linewidth 
of typically a few MHz. 


III. CRYSTALLIZATION 


While the photonic BEC regime has been explored 
previously for non-interacting photons thermalized via a 
dye [38] and in exciton-polariton condensates |39| , many- 
photon systems with strong interactions are wholly un¬ 
explored. Unlike cold atoms, for which interactions over¬ 
whelm kinetic energy only in the presence of an optical 
lattice to increase the effective mass of the atoms, the 
ability to tune the photon mass via the degeneracy of 
the resonator means that for realistic parameters, the Ry¬ 
dberg interactions can overwhelm the kinetic dynamics 
of the polaritons, allowing for studies of emergent crys¬ 
tallinity. In a near-planar cavity, the mass and harmonic 
oscillator frequency have a fixed relation to each other. 
However, for more general near-degenerate cavities, it 
is possible to design the harmonic oscillator frequency 
and the mass independently through the choice of cavity 
waist. 

Crystallization is expected beyond a critical value of 
the ratio between the interaction energy and kinetic 
energy in the system m- To obtain conditions for crys¬ 
tallization, we assume that the separation d between 


dark polaritons exceeds the blockade radius Then the 
effective interaction simplifies to Ud ~ —ocC^jp^ and 
Td = a\C&\mi:)/{h?d‘^). Crystallization therefore occurs 
at high density, as in the case of dipolar interactions. 
The density depends on the strength of the driving laser 
and the harmonic confinement. For self-consistency, and 
to prevent losses due to detuning from the EIT condition, 
the critical inter-particle spacing must exceed 


_ / amD\^ C&/h 
rdh ) \hx\^ 


( 9 ) 


The relation gives a condition on the polariton mass 
and therefore on the degree of cavity degeneracy for 
a given cavity waist. To obtain concrete values, we need 
to know Td at the critical point. Although this quan¬ 
tity is currently unknown for van der Waals interactions, 
quantum Monte Carlo simulations m of crystallization 
in two dimensions with dipolar interactions give a critical 
value of Td = 18 ± 4; we use = 20 as an estimate here. 
For the parameters given in Fig. [3] using the 80S Rydberg 
level and A = 0, the condition H then gives wj_ < 5.8 
MHz, which current technology for optical cavities can 
readily achieve. 

The emergent few-body structures for a nearly degen¬ 
erate cavity with s = 2, allowing only even harmonic os¬ 
cillator wavefunctions, are shown in Fig. [^ The Hilbert 
space restriction for s = 2 results in a unique type of crys¬ 
talline order, as illustrated in Fig. Eh- The symmetry of 
the photonic modes implies that photons localized away 
from p — 0 occupy superpositions of two locations, and 
these superpositions must fit in with the crystal struc¬ 
ture. As a basic check of this idea, we numerically sim¬ 
ulate the driven cavity system for up to two photons, 
and find that when the interaction strength is sufficient, 
the photons become localized to different positions, in- 
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FIG. 5. Laughlin states of light, a, Probability 
for two emitted photons to occur in the two-boson Laugh¬ 
lin state, versus the collective cooperativity ?7coii = /{FT), 

with photons injected into the I = 1 mode, b, Decompo¬ 
sition of the two-photon wavefunction into the basis 1 / 1 ^ 2 ) 
with the angular momentum indices of the two pho¬ 
tons. The states |11) and |02) are replaced with the states 
\L), \A) — T |02)), giving the Laughlin (-) and anti- 

Laughlin (-f) states, respectively. The color scale gives the 
probability for two photons to be in the given state. 


dicating the two-photon analogue of crystallization. The 
three peaks indicated in Fig. [^-c result from one photon 
localizing at p = 0 with the other localizing at two sym¬ 
metric positions. Figure]^ shows the ground state of a 
two-polariton system while Fig. & shows the normalized 
two-photon part of the steady-state wavefunction for the 
driven system. The calculation in Fig. & uses the same 
parameters as given in the caption to Fig. [^except that 
in 4c, Vl = 2 tt ■ 2.5 MHz and uj±_ = 27r • 500 kHz. For clar¬ 
ity, Fig. & uses a larger interaction strength = 0.7wo) 
and weaker harmonic trapping w_l = 27r-5 kHz to increase 
the separation between photons. 


IV. LAUGHLIN DROPLETS 


For a cavity with Landau-level-like degenerate man¬ 
ifolds, interactions between dark polaritons lead to 
bosonic Laughlin states [33 |S3| as many-body eigen- 

modes of the system. Unlike solid-state realizations |S3] 
and cold-atom proposals m 156] . where the Laughlin 
state is realized through thermal phase transitions and 
adiabatic state transformations, respectively, the open¬ 
ness of the photonic system enables spectroscopically- 
resolved preparation of the photonic Laughlin state. 

Laughlin states arise in the limit of a small blockade 
radius ^ <C ruo where the effective interaction reduces 
to a contact interaction. Additionally, we consider the 
limit of small mode splitting a;_L —>■ 0 to reach a Landau 
level. When photons are injected into the lowest Lan¬ 
dau level with angular momentum I — 0, strong blockade 
allows only a single polariton into the system, owing to 
the repulsion between the intracavity polaritons. When 


photons are injected into the I = 1 mode of the lowest 
Landau level, the blockade is lifted for two photons, as 
they are able to enter the system in a bosonic Laughlin 
state [53] . 


\^2l) — 




^ A'^ 

^2 ^0D^2D 


| 0 ) ( 10 ) 


with the creation operator for a dark polariton in 

mode m. In position space, (101 has the spatial wave- 
function 


(Pi, P 2 IVI/) OC (zi - (11) 


with 12 = 1/2. The Laughlin state has vanishing contact 
interaction energy, and non-zero overlap with two pho¬ 
tons in the I = 1 mode. Figure shows results from a 
simulation of the driven system in the two-photon Laugh¬ 
lin state setup. The overlap of the normalized two-photon 
part of the wa vefu nction \^ 2 ph) with the two-photon 
Laughlin state (10) is given by = | {'^ 2 ph\ ^ 2 l)P- By 
varying the atomic density up to 1 pm“^, we see that 
as the collective light-matter coupling is increased, the 
state of two photons in the system approaches the Laugh¬ 
lin state with high probability. Coupling photons instead 
into the I = N—1 mode of the lowest Landau level would 
lead to creation of the iV-boson Laughlin state with total 
angular momentum L = N(N — 1), due to the vanishing 
blockade in that state. In the Appendix we generalize 
the description of photonic Laughlin states to manifolds 
with s > 1. 

In conclusion, we have demonstrated that Rydberg po¬ 
laritons in multimode optical cavities provide a powerful 
platform for quantum simulation of strongly interacting 
photonic systems. The dark polaritons are well-described 
as particles in two dimensions and their Hamiltonian 
can be controlled through the design of the optical cav¬ 
ity. In particular, degenerate manifolds of Hermite-Gauss 
modes naturally lead to a harmonic oscillator Hamilto¬ 
nian, while degenerate Laguerre-Gauss modes in helical 
cavities naturally give rise to an effective gauge field. 
The interactions between dark polaritons can be strong 
enough to induce single-photon nonlinearities, allowing 
simulation of interacting quantum systems beyond the 
mean-held regime. In particular, we have shown that 
this platform allows studies of emergent crystallinity and 
fractional quantum Hall states. Our work motivates ex¬ 
periments on this system, as well as future theoretical 
studies of the system. In particular, this system gives ac¬ 
cess to many-body physics in a driven, dissipative regime, 
and promises to extend our understanding beyond con¬ 
ventional equilibrium and near-equilibrium studies. 
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Appendix A: Theory of Degenerate Cavities 

The frequencies of cavity modes, within the paraxial 
approximation, are given by [32103] 


waist, di is the detuning of the probe laser from EIT 
resonance, and ij^a is the photon field operator, defined 
in greater detail below. We work in the limit of weak 
driving, i? —>■ 0. In this limit, the total probability of an 
excitation in the cavity becomes vanishingly small, mak¬ 
ing the average time between quantum jumps arbitrarily 
large. The steady state of the system then approaches a 
pure state m- 

We obtain the steady state wavefunction |4's) by per¬ 
turbation theory in the weak drive U, 


= — [{m+ l/2)xi + {n + 1/2)x 2 + 27rp] (Al) 


where m and n are positive integer transverse mode in¬ 
dices, p is the integer longitudinal mode index, c is the 
speed of light, and Lr is the round-trip cavity length. 
The Gouy phases xi ,2 are given by the four eigenvalues 
g±*xi ,2 Qf round-trip ray matrix. 

Degeneracies occur when the Guoy phases are tuned 
to rational fractions of 27r. To obtain a two-dimensional 
harmonic oscillator Hamiltonian, one uses Xi = X 2 = 
27rq/s -|- e. The mode frequencies (Al) within a nearly 


degenerate manifold are then given by 


with 


« (l + GoU + GoUGoU)|0) (B2) 


Go{Sl) = 


- H 


Q, 


(B3) 


where |0) is the ground state of the atom-cavity system, 
with no photons in the cavity and all atoms in the ground 
state, and Q = 1 — |0) (0| projects onto the space orthog¬ 
onal to |0). Here H is the atom-cavity system Hamil¬ 
tonian not including V, and N^xc counts the number of 
excitation quanta in the system. 


= w_L(m + n -I- 1) -I- const 


(A2) 


with u}± = ccjLr-, and the constraint that m -I- n mod s 
is constant within the manifold. 


For a helical cavity, the indices m and n in (Al) count 


units of positive and negative orbital angular momentum 
in the transverse plane. Landau levels are realized when 
one of the Gouy phases, say xii is tuned to a rational 
fraction of 27r. The case Xi 0 results in conventional 
Landau levels containing a state at each value of the or¬ 
bital angular momentum. In general, for xi = 27rg/s, 
each degenerate manifold includes modes with orbital an¬ 
gular momenta in steps of s. Tuning slightly away from 
degeneracy, so that xi = 2'Kqf s + e, induces a harmonic 
potential in addition to the effective magnetic field. 


Appendix B: Numerical Simulation 

To simulate the two-photon crystal and Laughlin 
states, we obtain the steady-state wavefunction of the 
system in the limit of a weak probing beam. The system 
is modeled using 30 atoms positioned randomly within 
the atomic cloud, which is sufficient to approximate the 
continuum limit. The coupling strength is scaled up to 
the appropriate value to simulate a given atomic density. 

We carry out the calculation using perturbation theory 
to second order in the weak probe field to capture one- 
and two-photon effects. The probe laser acts through a 
Hamiltonian term 


Nexc = X! j (B4) 


The states resulting from (B21 have up to two excita¬ 


tions. Note that the apparent pattern in (B2) is modified 


beyond second order by the non-zero second-order shift 
of the eigenvalues in a discrete system. 


Appendix C: Laughlin States in Generalized Landau 
Levels 


The analog of the lowest Landau level (LLL) consists of 
the manifold of modes with angular momentum quantum 
numbers 0, s, 2s, etc. for some positive integer s. The 
conventional LLL corresponds to the case s = 1. 

We first consider eigenstates of the low-frequency ef¬ 
fective Hamiltonian for two dark polaritons with total 
angular momentum 2s. The wavefunction is of the form 


| 4 -) = 


(^sd) + b^0D^2s,D 


| 0 ) 


The eigenstate with zero interaction energy has 
a/b = —s!\/2/(2s)! 

In position space, the wavefunction becomes 


(Cl) 


(G2) 


(G3) 


V = HR j dp [vL{p)e +h.c.] (Bl) 

where R is the driving strength, fz,(p) is the mode profile 
of the probe laser in the transverse plane of the cavity 


up to normalization, where Pj = {xj,yj) is the position of 
the j-th dark polariton and Zj = xj -I- ipj. This state gen¬ 
eralizes the two-boson ly = 1/2 Laughlin state. To gener¬ 
ate this state experimentally, one can couple photons into 










the cavity mode with angular momentum s; interactions 
then populate the modes with angular momentum 0 and 
2s. Two-polariton states with total angular momentum 
ns for even n can be realized by coupling into the mode 
with angular momentum ns 12, resulting in the analog of 
the n = 1/n two-boson Laughlin state oc {z{ — z|)". 

For N dark polaritons with total angular momentum 
L = N{N — l)ns/2, the eigenstate in position space is 

{p„p„ ...,p^\^) cx n (C4) 

i<j k 

up to normalization. To realize this state one must couple 
into the cavity mode with angular momentum L/N = 
(iV- l)ns/2. 


which acts as a delta function on the space of functions 
spanned by \vmn}(m.n)&M- 

The field operator ( |D3| ) allows us to write the second- 
quantized Hamiltonian for photons in a given manifold, 


Hph — ^ ^ hitj- 


mn^rnn^'^n 




= j dpipl hph tpa 


(D5) 

(D6) 


using (D2| and (D3|, and defining the first-quantized 
photon Hamiltonian using hphVmn = ^mnVmn- 


Appendix E: Two-Dimensional Atomic Excitation 
Operators 


Appendix D: Cavity Photon Field Operator 

While a photon in an optical cavity moves in a three- 
dimensional space, the degrees of freedom become two 
dimensional after restricting to a nearly-degenerate man¬ 
ifold M of the types described in the main text. Two di¬ 
mensional behavior arises because the longitudinal mode 
index p becomes a function of the transverse mode in¬ 
dices m and n within the manifold. The photon then has 
transverse degrees of freedom but no longitudinal degree 
of freedom. 

To define the photon field operator, we first define pre¬ 
cisely the transverse mode functions Vmn(x,y). In gen¬ 
eral, the cavity modes are described by orthonormal elec¬ 
tric field amplitude functions Umnp{x, y, z). Near a waist, 
the mode functions factorize into a transverse part and 
a longitudinal part. 


In this section we provide a precise definition of the 
two-dimensional collective operators for the atomic exci¬ 
tations. To start, we consider the field operator 'I'„(r) for 
an atom at position r in internal state a. As the atom- 
photon interactions do not create or destroy atoms, but 
simply change their internal state, we define the local 
excitation operators 

= '^Ux)'^gix)/\Jng{r) (El) 

that promote an atom at position r from the ground state 
to the state a. Here ng{r) is the density of ground state 
atoms at position r. We work in the limit of a large 
ground state density and a small excited state density, 
so that ng, and 'I'f become classical numbers and 
d/J —>■ In this limit, the excitation operators (pa 

obey the commutation rule 


Umnp{x,y,z) Vmnp{x,y)Wmnp{z) 


(Dl) 


which remains valid within the Rayleigh range of the 
waist. Within a given manifold M, as p depends on m 
and n, we write the transverse part simply as Vmn(x,y). 
The functions Vran{x,y) are chosen to be orthonormal, 

J dp dp V,^„^{p)Vm'n'{,P ) — ^mm' ^nn' (D2) 

The two-dimensional photon field operator restricted to 
the manifold M is then defined as 

V’a(p) = E ^mn ip) Ocmn (D3) 

where amn is the annihilation operator for a photon in the 
(to, n) mode in manifold M. The ladder operators have 


1 


— ^m m.' 


mm '-'nn 


the usual commutation rule 
The commutator of the photon field operators is then 

[^a(p),V’l(p')] = E ^rnniP)v^rnnip') (D4) 

{m,n)£M 


(/)a(r),((ij,,(r') = (5„a'(5(r - r') 


(E2) 


To define collective excitation operators, we look at the 
atom-photon electric dipole interaction, 

Hed=^ f dr ^ ^ng{z)(j)l{r)umnir)amn+i^-c. 

(E3) 

where Iigo/2 = de\/ with Wg the resonant frequency 


2eo 

of the ground to excited state transition, dg the dipole 
matrix element, and cq the electric constant. The atomic 
density is assumed to be uniform in the transverse dimen¬ 
sions so that rig = ng{z). Using the factorization (Dll, 
Hgd can be written in terms of the two-dimensional pho¬ 
ton field operators (|D3|) as 


Hed = ^ J dp dz ng{z)4>\{p, z)w{z)tpaip) + h.c. 

(E4) 

where we have used the condition that the atomic sample 
is thin compared to the Rayleigh range so that the factor¬ 
ization (Dl) applies over the relevant range of integration 













9 


in z. Additionally, we have introduced the assumption 
that the longitudinal part of the cavity mode is indepen¬ 
dent of the mode index, so that Wmn{z) = w{z). This 
condition holds in a near-planar two-mirror cavity and 
in travelling-wave cavities with any value of s, but not in 
two-mirror cavities with degeneracies having s > 1, such 
as confocal cavities. 

After introducing the two-dimensional field operator 
for atomic excitations. 


looking at the Hamiltonian term for the coupling from 
the intermediate state to the Rydberg state in the rotat¬ 
ing wave approximation, 

H^ = h^ J dr + h.c. (E12) 

where kc is the wavevector of the control beam. We de¬ 
fine the two-dimensional field operators for the Rydberg 
level as 


'4’tip) = J \jz) (E5) — J dz Wj{z)e^^'^z) (E13) 


the electric dipole Hamiltonian (E4| further simplifies to control field coupling ( E12 ) becomes 


Hed = ^ Vh J dp tfjl{p)i;a{p) + d.c. (E6) 

The requirement that ■i/'e satisfies the commutation rule 
[^e{p),^l{p')] =S{p-p') (E7) 

fixes the value of the constant n at 

h = y dz ng{z)\w{z)\'^ (E8) 

For a sample with thickness d in the 2 direction, h ^ 


Ugd/L. Equation (E6) shows that the collective Rabi 


frequency is given by G = go'/^. 

While cavity photons in a given manifold only excite 
atoms into collective states with a particular longitudi¬ 
nal mode structure, the atomic gas supports arbitrary 
collective excitations. Additional atomic field operators 
account for these modes. Defining wq = Wy/ug/n, con¬ 
sider a complete orthonormal set of functions {u)j}'jLo 
and let 

^Ijip) = / Wj{z)(l)lip, z) (E9) 


so that V’eo = V'e from (E5). The 'i/'ej operators allow 


expression of arbitrary three-dimensional operators. For 
example, the density of atoms in the intermediate state 
is expressed as 




(E14) 


3 = 0 ' 


The two-dimensional field operators ipf.j and iprj satisfy 
the commutation rule 


= Sca'Sjf6{p - p') 


(E15) 


with a, a' = e, r. 

The Hamiltonian for Rydberg EIT in a nearly- 
degenerate multimode cavity can then be written in 
therms of the two-dimensional excitation operators as 


H = h Ip, 


^a\ ^ ( V G/2 0 



G/2 A D/2 
V 0 D/2 -i'y/2 

f A D/2 
Iprj J \ D/2 -17/2 



(E16) 


with hph = hph/fi — inj^. The second line in (E16) 
accounts for matter modes with longitudinal structures 
that do not couple to the cavity. These modes contribute 
to the screening of the polariton-polariton interaction. 


Appendix F: Diagonalization of the Single Polariton 
Hamiltonian 


J dr Pipe = J dppljpJej (ElO) 


3=0 


where we have use the completeness relation 


3=0 


(Ell) 


The single-particle eigenstates of the system Hamilto¬ 
nian (E16) can be found by working in a basis defined by 


the cavity modes. To express atomic excitations that do 
not couple to any cavity mode in the chosen manifold M, 
we extend the basis defined by M to a complete orthonor¬ 
mal set of modes. The collective operators for excitations 
coupled to an arbitrary mode m are then defined as 


The EIT control field couples the collective excitation 
generated by absorbing a cavity photon to a specific col¬ 
lective excitation of the Rydberg state. We therefore 
define two-dimensional collective operators 'pr{p) for Ry¬ 
dberg excitations in a manner analogous to that of the 
intermediate state. We identify the correct operators by 


Dm ^ J (t^eUm (FI) 

for the intermediate state, and 

Fm= J dr (F2) 
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for the Rydberg state. Here «„(!■) = Um{^)y/ng/n for 
m S M. For m ^ M, the Um are a basis for the or¬ 
thogonal complement of the cavity modes in the space 
of normalizable functions, so that {um} taken over all m 
form a complete, orthonormal basis. The collective mode 
operators obey [Dm,Dl] = [Fjn,FX\ = S^n- 

The fundamental form of the system Hamiltonian is 
given by 



Fffint (F3) 


where Sm is the detuning of mode m from EIT resonance. 
Using the collective operators (FI) and (F2|, the non¬ 
interacting part of the system Hamiltonian Fq = FI — 
Flint can be written as 


a. Dark Polaritons 


We diagonalize the 3x3 matrix in (F4) as = 

QmEmQTTi 5 where The 

columns of are right eigenvectors of H,„ while the 
rows of are left eigenvectors of H^. The annihi¬ 
lation operators for dark and bright polaritons in mode 
m are then 4l,72_|_, — Qm with 

Jm = (cjn, Dni, Fm)'■ The dual creation operators are 

= JLQm- The polariton oper¬ 
ators obey the commutation rules = 5mn5p.v 

and = 0 . 

Field operators for the polariton excitations are 
then defined as = Y^maM'^rnAmii and = 

^r fi € {D,+,-}. The commu¬ 
tator for the field operators is ['ipf_i{p),tl}'^{p')] = 

the sum acts as a delta func¬ 
tion on the space of functions spanned by {vm}- 

The term in the Hamiltonian (F4) describing the dark 


polariton excitations then becomes 


Ho = 




G/2 0 \ / 

A H/2 

n/2 -fo/2 / \ 

n/2 \ {Dm\ 

-i7/2 J \FraJ 


dm \ 

Dm 

Fmj 


(F4) 


Hd = ^ AfnD^mDAmD (F 6 ) 

mGM 

= J dp tj}f)hdark{p)lljD (F7) 

where /idark(p) is a differential operator on functions of 
the transverse coordinates p = (x,y) such that 


The single-excitation eigenstates of H in the photonic 
modes m G M are familiar EIT eigenstate s I46j . ob¬ 
tained by diagonalizing the 3x3 matrix in (F4). The 


two bright states have complex eigen-frequencies SmB± 

f A ± \/A2 + G^ + 


/ 2 , to zeroth order in Sm and 7 . 

The dark state eigen-frequencies SmO = dmo — i^mol^ 
are given by 


OmD 




G 2 


{K-i'f) 


4G2H2(A-i0 

(G 2 + H2)3 


G' 


G 2 

.7 — 


(--i) 


(F5) 


to second order in the small quantities Sm, Km, and 7 . 
The leading order probability for a dark state polariton 
to be a photon is cos^ 6 = n^/{G^ -I- where 9 is 
the dark state rotation angle. The first term in (F5) 


therefore gives the probability for the dark polariton to 
be a photon, times the complex detuning of the photon, 
while the second term gives the probability for the dark 
state polariton to be a Rydberg excitation, tim es the loss 
in the Rydberg state. The second line in (F5| gives the 


second order eigenvalue, reflecting modifications to the 
dark state wavefunction due to detuning and losses. In 
particular, a non-zero detuning Sm from EIT resonance 
leads to population of the atomic intermediate state and 
to a loss rate proportional to TS^. 


^dark^m — ^SmD’^m 


(F 8 ) 


From (F5) we see that /idark is given to leading order 
by ^dark = feph cos ^ 9, with kph the photon Hamiltonian 
introduced in (|D 6 |). For the harmonic oscillator-like and 


Landau level-like manifolds of cavity modes. 


^dark(p) — 


1 


1 


2m D 


{-ihdp - A) -h , (F9) 


where the effective gauge field A ^ 0 for the Landau 
level-like case. Note that the polariton mass is increased 
relative to the photon mass by the factor 1 / cos^ 9 and 
the osci llati on frequency is decreased by cos^ 9 so each 
term in (F9) equals the corresponding term in the photon 


Hamiltonian, scaled by cos^ 9. 


Appendix G: Two-Polariton Problem 

The effective interaction between dark polaritons fol¬ 
lows from the solution of the two-body problem 

= (iLo+i7i„t)|«') (Gl) 

where is an arbitrary two-excitation state. 

The problem can be re-written as 

14-) =Go(a;)iLint|^) (G2) 
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where 

Go{u^) = [nw-Ho]-^ (G3) 

The inverse exists for most values of uj because this is a 


with, 


C/®® = ^—^-TTT and (G13) 


i-xv 


1-x^ 


discrete system. Equation (F4) shows that Hq is block 
diagonal in pairs of modes (m, n). We find the Green 
function (G3) by inverting within each block. 


Equation (G12) implies Pip = 0. Therefore we expand, 


As the interaction term i/i„t projects |d>) onto the 
subs pace containing two Rydberg excitations, projecting 
(G21 onto the two-Rydberg subspace gives a Schrddinger 


equation for Rydberg-Rydberg part of the wavefunction 
alone. 


‘P(x,y) = ^ X! + 6'm(x)u™(y)] , (G14) 

m^M 

with 9m as arbitrary functions, 

0m{r) = ^amnUn{Y) (G15) 


p’rr — GrrP'mtp’r 


(G4) 


where iprr is the Rydberg-Rydberg component of and 
Grr is the component of Go that takes two Rydberg ex¬ 
citations to two Rydberg excitations. In the mode basis, 
Grr is diagonal in pairs of modes. Its elements consist of 
three types, corresponding to three cases for the pairs of 
modes: neither in M, one in M and one not in M, and 
both in M. We decompose Grr into these three types 
using projection operators, 

Grr =XP+ Y. XliSm)Pm + Y mj 

mGM m,nGM 

(G5) 

where Pmn projects onto the pair of modes m and n. In 


The Schrddinger equation (G12) can be rewritten as 


an algebraic equation for the expansion coefficients amn, 
giving 

aki = [xi{5k)-x] Y (G16) 

m£M;n 

for k G M, I ^ M and 

aki+aik = [xo{Sk,6i)-x] Y + 

m^M\n 

(G17) 

for k, I G M. The matrix elements of the effective poten¬ 
tial are defined as, 


the position basis. 




Pmn{x, y; x\ y') = Um{,x)Um{x')Un{y)u*r^{y') 

(G6) 

= J dxdy Mfc(x)u;(y)C/®®(x-y)'u^(x)M 

n{y) 

in terms of which the two other projectors are 


A nm 
— ^Ik 

(G18) 

Pm — ^ ^ {Pmn Pnm^ 

(G7) 

b. Low-Energy Limit 


n^M 



P — P 

J ^ mn 

(G8) 

When the interaction is sufficiently weak, or when 

m,n^M 


there are sufficiently many modes so that the 

excita- 

The projectors satisfy the completeness relation: 


tions can be spatially separate, the Hamiltonian has small 
eigenvalues oj. In the a; —)■ 0 limit, Xo{^m,^n) becomes 

P + Pmn = 1) 

mGM m,nGM 

(G9) 

large for small Sm, Sn and dominates over the other terms 
Xi{^k) and y. Equations (G16) and (G17) give 

where 1 is the identity operator on the space 

of two- 

(tkl 0 

(G19) 

Rydberg wavefunctions. 

We obtain the effective potential by subtracting yl 

for k G M, 1 ^ M and 


irom Grr? 


^kl P ^Ik ^ Xoi^ki ^ ^ i^mn P ^nm)^kl 

(G20) 

AG = Grr - xl 

(GIO) 

m^n£M 


= Y lXl{Sm)-X]Pm 

m£M 


Eor k,l G M. Also, ip becomes 


P ^ ^ [Xo(^m? X] Pmn 

(Gil) 

ip{x,y)-i'^ Y {amn + anm)Um{x)Uniy) 

(G2I) 


m,n^M 


The two-Rydberg component (G4) of the Schrddinger 
equation can then be re-written 


m,n^M 

So only the sum Umn + o-nm has significance. Define 


(^ = AG 


(G12) 


^mn — 2 “f ^nm) • 


(G22) 
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Also note that the asymptotic form of xo is 


c. Polariton-Polariton Interaction 


Xoi5k,Sl) 


with a = 


q4 


£k- 

(G23) 

and ek the eigenvalue of the single-polariton Hamiltonian 
for a dark polariton in mode k. 

The low-energy Schrodinger equation (G201 becomes 
w bki = (efc + £i)bki + O' 'y ' A^i'^bmn (G24) 

m,nGM 


The low-energy effective Schrodinger equation (G24) 
describe excitations in the cavity modes interacting via 
the effective potential 


Finally we show that the full eigenstate ip simply de¬ 
scribes a pair of dark polaritons in a superposition of 
cavity modes, via 


Pj = (G25) 


eff 


So far we have obtained an effective potential U' 
three dimensions. To define the effective potential be¬ 
tween dark polaritons we note that the interaction ma¬ 
trix elements of Eq. (G18) can be written, when all 
the indices refer to modes in M, as 


= j dpdp' vl{p)v*i{p')W^{p- p')vm{p)vn{p') 

(G26) 

with 

U‘^^{p) = J dzdz'\w{z)'w{z')\^ {p, z — z') (G27) 

In the limit of a sample that is thin compared to the 
blockade radius. 


^ [/®®(p,0) (G28) 

Finally, the interaction between dark polaritons is 

Ud{p) = aW^{p) (G29) 


which requires the use of all rows of the two-Rydberg col¬ 
umn of each block of Gq. In the limit where the Ck are 
smaller than the is a function of the restriction 

of C/®® to M. Since Gq suppresses the non-M compo¬ 
nents, as well as the non-dark polariton components, the 
end result of (G25) is a superposition of pairs of dark 
polaritons in the same modes described by (p. However, 
when the are larger, t/®® mixes the modes of p. 


We can then write the interaction between dark polari¬ 
tons in second quantized notation. 


Hd 


J' (G30) 

+ ^ / dpdp'p^l){p)pjl)ip')UDip - p')^Pd{p)^d{p') 
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